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1. Introduction 

After the considerable success of two dimensional quantum field theory and statistical 
mechanics, - integrable models on 2D regular lattices, conformal held theories, Liouville 
theory and matrix models of 2D gravity and non-critical strings - progress in analytical 
results in this held has slowed down. 

Among the principal questions remaining unsolved are, hrst, the so called c = 1 
barrier for non-critical strings (c is the central charge of the matter), and, second, the 
mysterious connection between the physical properties of various integrable 2D models 
coupled and non-coupled to gravity. The hrst problem is usually atributed to the absence 
of a stable vacuum for c > 1, though it has never been clearly understood. Indeed, in 
terms of matrix models, the obstacles seem to be purely technical. The second problem 
concerns the observation of many intriguing relations between 2D physical systems with 
and without coupling to 2D quantum gravity. 

This relation is clearly established on the level of critical exponents by the use of 
the continuous formulation of 2D gravity [|l[. The conformal dimensions of matter helds 
undergo a simple quadratic transformation as the result of gravitational dressing. The 
same phenomenon is observed, of course, in the matrix model formalism. However, if 
one goes away from the critical point, the relation between the physical properties with 
and without gravitational coupling, although persisting, becomes much more tricky and 
fragmentary. For example there is the description of 2D gravitational systems in terms 
of KDV hierarchies of classical 2D integrable systems [|| , as well as a strange coincidence 
between the amplitude for open string in SOS formalism and the S-matrix of the two 
dimensional Sine-Gordon theory [^. 

It seems that two-dimensional physics is more united than one would think at hrst 
sight. An interesting question to ask would be the existence of some interpolating models 
connecting the gravitational and “hat” phases of the same matter helds. 

Our paper is inspired by this physical idea, though it concerns mostly the elaboration 
of a technique for the solution of a new type of matrix model. The model describes, in the 
large N limit, planar graphs having arbitrary coordination number dependent weights for 
both the vertices and faces. In other words, we introduce a set of couplings ti, t2, ■■■tq ,..., 
the weights of vertices with 1, 2,..., q ,... neighbours, and t^, ..., the weights of the 

dual vertices (or faces) with appropriate coordination numbers. In Fig. 1 below is a 
typical graph with, for illustration, some vertices on the original and dual lattice labeled 
with their associated weights. The matrix models under consideration allow us to generate 
the following partition f un ction of closed planar graphs G: 

n (I'l) 

G Vq,V*eG 
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where Vq,v* are the vertices with q neighbours on the original and dual graph, respectively, 
and ^v* are the numbers of such vertices in the given graph G. We propose to call 
this the model of dually weighted graphs (DWG). 



Fig. 1: A typical surface and some of the associated weights 

It is clear that this model opens the way to understanding the very interesting transi¬ 
tion mentioned above. If we set ^4 = = 1 and tq = t* = 0, for Q 7 ^ 4, only regular square 

lattices (graphs) will exist in equation (| 1 . 1 |) . Hence, there are trajectories in the coupling 
space of this model, interpolating between pure gravity (for example, when all t* = 1 ) and 
the regular “flat” lattice. 

We will show in this paper, that the underlying matrix model describing the DWG is 
solvable. Our solution is based on an elegant representation of this model in terms of the 
group character expansion found in |^ . It allows us to reduce the N‘^ degrees of freedom of 
the original matrix model to the N degrees of freedom labeling a representation. We then 
apply the saddle point approximation to hud the most probable group representation in 
the corresponding sum over characters, specihed by the distribution of its highest weights. 
A similar approximation was hrst successfully used in [0 for the calculation of the QCD 2 
partition function on the sphere. We conclude with a well dehned (though complicated) 
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integral equation for this distribution. Though we have not yet been able to extract the 
physical picture corresponding to the “flattening transition” we demonstrate on a number 
of simpler examples that our method is consistent and correct. 

Furthermore, we solve a model apparently inaccessible by standard methods: we 
calculate the number of planar graphs having only even number of neighbours for both 
original and dual vertices t 2 q = = 1, and t 2 q-i — *?■ 

We hope that our methods will lead to new progress in solving many physically inter¬ 
esting 2D systems. A natural step forward would be the introduction of matter on DWG, 
a tempting opportunity, whose success is, of course, not automatically guarantied. Since 
our matrix model is a generalized matrix external held problem, it could also be useful for 
new studies in random (mesoscopic) systems. 

We present below explicit details of our technique as we feel it is a general and powerful 
method for matrix models. 


2. Reduction of the DWG model to a sum over characters 

The partition function for the dually weighted graphs can be formulated as the fol¬ 
lowing matrix model (see for example @ ): 

Z((,r) = A-^ j VM p J) 


The matrices A and B are external matrices. In the perturbative expansion of the above 
integral, the matrix B weights a vertex of coordination number q with the factor Tr B'^ 
while each face bounded by q vertices is weighted by a factor Tr We can therefore 
make the connection 

t, = i i TA and (J = i 1 TV A". (2.2) 

q N ^ q N 

Note that it is impossible to solve the above matrix integral by standard methods since it 
is unclear, for A ^ 1, how to perform the angular integration, ie. how to evaluate unitary 
matrix integrals of the form f (dO)//exp(^^/?fcTr(MfiAfil)^) with M and A diagonal. 
We circumvent this difficulty by expanding the potential in terms of the characters on the 
group: 


Tr{MAA ^ TT ^ 

(1 - R.(MA),) 


XRiMA). 

R 


(2.3) 
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Here Bi and {MA)j are the eigenvalues of the matrices B and MA. The hrst step involves 
rewriting the sum over /c as a double sum over all the eigenvalues of the matrices B and MA 
of — ln(l — Bi{MA)j). Exponentiating the log then gives the product in the numerator. 
The second step uses a group theoretic result to write the product in terms of a sum over 
characters. The character is dehned by the Weyl formula: 


X{h}{A) 


A(a) 


(2.4) 


where the set of {h} are a set of ordered, increasing, non-negative integers, A(a) is the 
Vandermonde determinant, and the sum over R is the sum over all such sets. The RA label 
representations of the group U{N) and the sets of integers, {h}, have the correspondence 
with the Young tableaux shown in Fig. 2. 



Fig. 2: Connection between young tableaux and the integers hi 


Note that the restriction on the allowed Young tableaux that any row must have at least 
as many boxes as the row below implies that the {hi} are a set of increasing integers 

/li+i > hi. (2.5) 

Substituting equation (|2.3D into the integral in equation ( p.l|) we can now do the angular 
integration using the identity f (T>0)j/Xi^(0M0lH) = xr{M)xr{A)/ dn (where is the 
dimension of the representation given by du = A(h)/ arrive at the expression 
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The gaussian integral can be done explicitly and we arrive at the hnal formnla 


Z = c 






■X{i)('4)X{i)(B)(-77) 


N- 


(2.7) 


where c is some nnmerical constant that we can drop, the {/i^} are a set of A^/2 even, 
increasing, non-negative integers and the {h°} are N/2 odd, increasing, positive integers, 
and the snm is over all snch sets. In other words the original snm is now restricted to the 
subsets of {h} with equal numbers of even and odd integers. This is an exact result. The 
sum is in general divergent, as is the original matrix integral, and should be thought of as 
a generating function for graphs of arbitrary genus. This formula was originally derived 
by Itzykson and Di Francesco by summing up “fatgraphs”. For the rest of this paper 
we restrict our attention solely to the genus 0 contribution. In other words we will be 


studying the large N limit of equation (|2.7|) . 

There is a second useful formula for the character given in terms of Schur polynomials, 
Pn{0), dehned by 

OO 

~ ( 2 , 8 ) 


n=0 


In terms of Schur polynomials the character is 


X{h}{A) = (2.9) 

where 9i = iTrZ*. In general, the explicit expressions for the characters are very compli¬ 
cated. Certain specihc cases however yield simple results which we will state as we need 
them. 


3. Relations between highest weight and matrix model quantities 

Before we look for the large N limit of equation (|2.7|) , we derive some explicit ex¬ 
pressions relating useful quantities from matrix models to quantities encountered in the 
language of highest weights. In the large N limit of equation (|2.7D , we assume that the 
sum over all representations will be dominated by a single contribution, or a single Young 
tableau, [hi], and introduce a density p{h) dehned in the standard way by p{h) = 

We are being slightly sloppy with the notation here since to dehne a sensible density we 
have to rescale the integers {hi} by dividing them by N. For the rest of this paper an hi 
with an index refers to one of the original integers and h without a subscript refers to the 
rescaled continuous parameter. 
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All the formulae in this section have their root in the simple observation that 


TrA-^ = 

k 


X{h}{A) 


where hi = hi + qSi^k- 


(3.1) 


For compactness of notation we have omitted labeling the h with an index k. This formula 
follows directly from the Weyl formula for the character, (|2.4|) . The character can be 
written in terms of the Itzykson Zuber integral, /(h, a) = det^i^ (e^''“‘)/(A(/i)A(Q!)), as 


Xh{A) = I{h, a)/X{h) 


A(a)’ 


(3.2) 


where ai is dehned through the eigenvalues of A by a; = . This allows us to write 


TrA'? = 


^(h) In J(^,Q;)-ln I(h,a) 

- -e^ 9 

A(ft) 


A” s n ( 1 + 


k j(^k) 


hk hj 


^qF{hk) 


(3.3) 


where 

F{hk) = j^logI{h,a) (3.4) 

ohk 

is the derivative of the Itzykson Zuber integral. In the last step we assume that the Itzykson 
Zuber integral has a well dehned large N limit: I{h,a) — F[p{h),p{oL)\+o{N°)_ 

notice that we can replace the sum by a contour integral in h (encircling all the hfc’s) if we 
also unrestrict the product allowing it to be a product over all j. This contour integration 
trick was originally invented in P] for a much simpler model (c = —2 gravity), and more 
recently was used in 0 for the large N limit of the heat kernel. Applying it in this case 
we obtain: 

— TrA-? = - / ^eF^(h)+F{h)) _ f (3.5) 

N q J 27ri J h — h' 


F{h) is initially dehned only on the support of h; we have then analytically continued F{h) 
into the whole complex plane, so that the contour integral, which circles the cut of the 
resolvent, H{h), is well dehned. We have thus rederived in a very compact way the large 
N limit of the Itzykson Zuber integral |^] [1^] . To make the co nn ection with the result 
in IP more explicit it is simple to expand (|3.5|) as a power series in q and then resum the 
series to express the result in terms of the resolvent for the eigenvalues a (see also ): 


Q{a) = — <j) — H{h) — F{h)) with ©(a) = J 


da'a{a') 


a 


a' 


(3.6) 


where a{a) is the density for the eigenvalues a. 
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Next we look at the expectation value of Placing this term into the M 

integrand in equation (|2.6|) and using similar steps as for ■^TrA’^ we arrive at 

—TrM^^ = — / (3.7) 

N q J 2Tii 

This formula is derived only for traces of even powers of the matrix M. The expectation 
value of ^Tr(MA)^'^ is derived in a like ma nn er. This time we substitute ^Tr(M74)^'^ 
into an earlier step in the derivation of the Itzykson-Di Francesco formula, namely into 
equation ( p.3|) . Using (|3.1|) we see that this time we shift both ^md xr{M), which 

leads to the hnal result 


^Tr((MA)29) = y y" 


Again this is derived only for traces of the even powers of MA. We assume from 
that we are working with an even potential so that it is only the even traces that 

over all q, assuming also that the potential 


Summing up equations (|3.7|) and 
we arrive at two formulae for the two types of resolvents one can dehne for the 
matrix model: 


here on 
remain, 
is even, 
original 


WiP) =(lTr 
WAiP) =(^Tr 


P-M^ 

- - - ) 

p-ma' 


1 

p ~ 
1 

“ p 


_ ^ ^ ln(p2 _ xhe^ih)+2Fih)^ ■ 


(3.9) 


The hrst equation can be solved by Lagrange inversion since we know that the only sin¬ 
gularity of H{h) is the cut circled by the contour. Performing the inversion gives the very 
simple pair of equations 


p2 

PW{P) =— - h 


(3.10) 


Knowing H{h) we perform a functional inversion to obtain the resolvent of the original 
matrix model. We cannot in general do the same inversion for the resolvent Wa{P) since 
we do not know in advance the singularities of F{h). 


4. The gaussian model: straightening random loops 


We will now check the power of our method on the simplest non-trivial case of our 
external held problem: We simply set B = 0, i.e. all tg = 0 in equation (|2.1|). Now there 
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is no potential at all and thus no weights are excited: The only contribution to equation 
( p77|) is the empty Young tableau: 




(4.1) 


One immediately checks that the inversion formulae (|3.10|) correctly give the Wigner semi¬ 
circle law (we may set A = 1 here) 


W{P) = \{P- \/P^ - 4). (4.2) 

Clearly W{P) and H{h) are unchanged even in the presence of non-trivial coupling con¬ 
stants t*, but now the interesting quantity is the resolvent Wa{P)'- 


Wa{P) = 1 


VM e 


-^TrM^ 



1 

P-MA' 


(4.3) 


We can hnd it by substituting H{h) into (p.5|) 


q = l 


= 


— In [ft - 1 - ft - 1. 

2tii ^ ^ 


(4.4) 


where we also summed up the moments constructing their generating function with an aux¬ 
iliary variable uj. One now sees that in this case the Lagrange inversion can be performed 
by picking up a pole term inside the contour, giving immediately 


OO 

h —1 = ^ qt* and h — 1 — u) h (4.5) 

q=l 


In addition, from the inversion formula (|3.9| ) for Wa{P) we hnd, using the same method, 
h = PWA{P) and h -1 = ^ (4.6) 

pz 

Eliminating F{h) we hnd the exact solution of our problem: 

CXD 

P^ QPq and P Wa{P) = 1 +P^ (4.7) 

q=l 


Indeed, this set of functional equations determines, for any set of couplings t*, after elimi¬ 
nation of u, the desired resolvent. Let us remark that these equations can alternatively be 
derived using Schwinger-Dyson techniques, yielding a non-trivial check of our functional 
methods. It is interesting to observe that we may obtain arbitrarily complicated resolvents 









by freely choosing the t*’s. On the other hand, it is seen that a hnite number of non¬ 
zero coupling constants always leads to an algebraic resolvent. An amusing toy system 
consists in only activating the hrst three coupling constants. Now the resolvent Wa{P) is 
interpreted as the generating f un ction of rainbow graphs with face-valency not larger than 
three, see Fig. 3 below. It is, from ( fl.7|) , given explicitly by 


Wa{P) = - + 


2t 


*2 


(P'^-t*f-2t*t;-J(P^-q)^-4tiq 


(4.8) 


It is interesting to investigate what happens in this toy system if we tune away the faces 
with negative boundary curvature, i.e. ts 0: 


Now we obtain merely the “cigar-like” diagrams below. 




Fig. 3: “Rainbow” —> “cigar-like” diagrams in the gaussian model 


In fact, as it was argued in Appendix A of [jT^, the universal continuum limit (with 
string susceptibility 7str = | due to the square-root singularity in equation ([4.8|)) of the 
model ( II.8|) can be interpreted as two-dimensional topological quantum gravity: The ex¬ 
pectation value of the metric tensor is zero in the bulk, leading to a theory of quantum 
loops. “Flattening” in such a theory is thus just “straightening”, and indeed the contin¬ 
uum limit of (|4.9|) , Fig. 2, is simply a straightened loop with two curvature defects. Note 
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that the cross-over from the “quantum” to the “straight” phase is simply a catastrophe in 
the algebraic sense: As soon as we tnrn on the negative cnrvatnre conpling tg, the defects 
proliferate and the straight line disorders. It is not exclnded that a similar rather trivial 
mechanism will govern the crossover from random to flat graphs. However, this is not the 
most likely scenario; indeed one is reminded of the fnndamental difference between one 
and two-dimensional systems with regard to the absence, respectively presence, of phase 
transitions. At any rate, let ns tnrn to real planar graphs. 


5. Saddlepoint equations and planar graphs 


Success was guaranteed in the case of the Gaussian model since we knew from the 
start the trivial, linear distribution of weights. We now have to establish that non-trivial 
distributions H{h) can be found and that we are able to reproduce planar graphs from the 
character expansion. One hnds that saddlepoint techniques may be successfully applied if 
certain precautions, to be elaborated below, are taken. Let us illustrate the method and 
its subtleties on a number of examples: 

In our new language, the simplest model generating planar graphs turns out to be the 
case A = 1 and B = J, where Tr equals one if q is even, and zero otherwise. Thus 
we obtain a traditional one matrix model with the “even-log” potential — | In (1 — M^), 
generating planar graphs with arbitrary even vertices. It is easy to explicitly work out - 
with the help of the Schur character formula ( p.9|) (see Appendix) - the characters for this 
case: 

X{/i}(l) ~ A(h) and X{h}{J) ~ A(h°) A(h^) sgn - h°). (5.1) 


From here on we will omit irrelevant numerical constants. One now sees that the character 
expansion (|2.7D for the partition f un ction becomes 


Z ~ A(h‘’)2 A(h^)2 (A) 2 (5.2) 


We thus observe that even and odd weights completely factorize! By symmetry, they should 
have the same statistical distribution. This partition sum is ideally suited for a saddlepoint 
analysis: the Vandermondes repel the weights from each other while the potential attracts 
- for small coupling A - to the origin. It is therefore natural to write down, in the large N 
limit, the saddlepoint equation 

/ 'i'*' >" ('* (5-3) 
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obtained in the standard fashion from equation (|5.2|) . The density p{h) and the continuous 
variables h were dehned in section 3 and one also uses Sterling’s formula: In h\\ ~ ■|(ln(h) — 
1). This equation is easily solved but leads to the wrong result. The phenomenon is 
identical to the one previously encountered in [0: The naive saddle point equation fails 
to take into account the constraint p{h) < 1 which follows from equation (|2.5| ). Imposing 
the condition that the density is saturated at its maximum value p{h) = 1 on the interval 
[0, b] , we write down the modihed saddlepoint equation 

f = (5.4) 


determining the non-trivial piece of the density on the interval [b, a]. We generate the full 
analytic function H{h) = f dh' p{h')/{h — h') from -f dh' p{h')/{h — h') by performing 
the contour integral 


H{h) 


ln(^^^) + ^{h-a){h-b) jf 


ds |ln(sA)+ln(^) 

27rz (s — h)-^(s — a){s — b) ’ 


(5.5) 


where the contour encircles the cut [a, 6]. Inflating the contour and catching instead the 
cuts [cx), 0] and [0, b] we arrive at 


H{h) 


In 


y/a — y/b h + a/^ + 

a/A {a-\-b)h — 2ab + [h — a){h — b) 


(5.6) 


with 




(5.7) 


the constants a and b being hxed by the condition that H{h) = l/h + 0{l/h‘^). This is the 
solution to the “even-log” matrix model in the language of highest weights, reproducing 
the correct critical coupling Ac = -1. Below we show a plot of p{h) for various values of A. 
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Fig. 4: Highest weight density p{h) for potential — ln(l — M^) 


It is easy to check the correctness of this solution by perturbation theory. Even 
better, we can independently calculate the eigenvalue resolvent W{P) of this model by 
traditional methods and demonstrate its exact coincidence with the resolvent obtained 
from the inversion formulae (|3.1C1|) . 

To demonstrate the power of our method we will next consider a case that is not 
a traditional one matrix model and thus has not yet been solved with other methods: 
Consider A = J and B = J, i.e. planar graphs with even coordination numbers for 
vertices and faces. Here eqs. (pi.7|) and (b-l|) lead to 


z~ 

i 






m) 


(5.8) 


Here even and odd weights no longer factorize. However, it is natural to assume that 
they are equally distributed. Thus the crossproduct should precisely cancel one power of 
a Vandermonde, leading to the saddlepoint equation 

This equation is solved exactly as the previous case and one hnds the weight resolvent 


H{h) = In 


(a — h) 2h‘^ 
hX{^/a + Vb)‘^ 


{y/a — a/6)^ + 2ab + 2{h + \fab) [h — a){h — h) 

{a + b)h — 2ab + 2v^\/ {h — a){h — b) 

(5.10) 
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with the interval boundaries being determined through the quantities ^ 

h=(^^)^by 

— ^ + 1 = 0 and rj =-{^ — 1) (5-11) 

O 

One easily hnds the critical coupling to be Ac = |. It is satisfying to observe that this is 
very slightly less than two times the value of the previous case; this is as expected since the 
asymptotic growth of the number of graphs with n edges is ~ A“"'. It is straightforward to 
verify that this solution indeed correctly co un ts the graphs under consideration; e.g. from 
equation ( b-lOP with the help of (|3.9|) : 


i TVM^ = A + A® + + 54A’' + .., 




1 
4 


-O' 






jfv 


A + 


A^ + O ^ A® + 




A^ 


+ ... 

(5.12) 


where the dots correspond to insertions of the matrix J'. 

These examples correspond to the ensembles of planar graphs simplest in the weight 
language. It is natural to ask for the description of the simplest original even model 0 


of pure gravity, i.e. the one matrix model with the action —^TtM^ + ^TrM^. It is again 
straightforward to explicitly calculate the characters here (see Appendix); the weights are 
now grouped into four blocks where e G {0,1, 2, 3} denotes their congruence, modulo 
four. This leads to the expansion 


Z-A-"^ A(hW)2 A(/i( 2))2 JJ(hf) 

{hO,h^} i,j 

^ A(h<‘))= A(h(«)^ n('T 

{h^,h^} i,j 


hf^) 4^an(A/Ai'V^)-l) 


(5.13) 


where for convenience we have substituted in Sterling’s formula for the factorials. One 
observes that even and odd weights factorize, but not the congruence classes (0, 2) and 
(1,3). In fact, each of the non-factorizing sectors has a structure identical to the case 
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A = B = 1 (i.e. the one matrix model with action —^TrM^ — ln(l — M)) since here the 


character expansion (|2.7|) gives, together with (|5.1|) , 


A-T ^ A(h”)2 A(h')" 


hA e 


— 


(ln(Xhk/N)-l) 


(5.14) 




Here it is the even and odd weights that remain conpled. We see directly that the partition 
fnnction in eqnation (|5.13|) is the square of the partition function in equation (|5.14|) . We 
thus rediscover in the highest weight language the well known connection between these two 
models. At the diagrammatic level this can be seen by placing a vertex of the model at 
the midpoint of every edge of the — ln(l — M^) diagrams so that the face centres of the 
model are the vertices and face centres of the — ln(l — M^) model. It is tempting to make 
a saddlepoint ansatz like in (|5.4|) , multiplying in this equation the principal part integral 
by an extra factor of | due to the variation of the extra factor ~ h^)- However, 

here the solution of this equation does not lead to the correct result. We can gain some 
insight into this failure by computing, by the usual means, the eigenvalue resolvent W (P) 
and then deducing H{h) from (|3.10| ). The result of this calculation leads to a third order 
algebraic equation for 


H(h) = ln(4h)) 


(5.15) 


with X [h) dehned through the solution of 


AA^ - A(1 + h)X^ + d-h + 7(A))A + h'^ = 0 (5.16) 

9 

with 7(A) = ^ (i-a/i- 12 A)(i- 12 A) ^ hnds on inspecting H{h) that the saddlepoint 

conhguration of weights is complex-. The rapid sign-changes of the product — h^) 

destabilize the reality of the saddlepoint. It is worth pointing out that the saddle point 
nevertheless exists, even though it is much harder to hnd. As we have seen in the previ¬ 
ous example, the presence of this factor in the denominator of the expansion is however 
without danger. A rough intuitive “explanation” is that in the numerator the product 
acts to repulse the different distributions, destabilizing the saddle point, whereas in the 
denominator it attracts and stabilizes. We will see in the next section that the stability 
of the saddlepoint can be preserved in the case of greatest physical interest: the gradual 
flattening of the random surface. 
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6. Flattening random lattices 


Before we flatten our surface it is worth understanding how the flat lattice is repre¬ 
sented in the language of highest weights. In this case ^^Tr = ^Tr and it is 

simple to derive the characters from equation (|2.9| ) (see Appendix), to obtain the partition 
function 


A“ 


E 


A(/iW)2A(h(i))2A(h(2))2A(/i(3))2 




( 1 ) 


hf){h\ 


( 3 ) uA)\(uA) 


hT){hr -hf){hf^ 




( 6 . 1 ) 


where again we have substituted in Sterling’s formula for the factorials. The potential term, 
gSfe,e 2 ^k attracts to the origin for A < 1, and repulses and is unstable for A > 1. The 

critical point is therefore now Ac = 1. The repulsion of the Vandermondes in the numerator 
is now precisely balanced by the attractive effect of the products in the denominator. 
Indeed, using our rule of thumb that the variation of a product Hi in the 

denominator precisely cancels the variation of one power of a Vandermonde, we arrive at 
the trivial saddle point equation 

-^ln(A)=0. (6.2) 

We observe that the saddle point equation is not satished anywhere (except possibly for 
Ac) i.e. none of the weights are excited and, as in the gaussian case, only the empty Young 
tableau contributes: 

«('*)= In (6.3) 

To order V^, only the original gaussian term can contribute to the partition function and 
expectation values of TrM'^'. The potential term cannot contribute to either, since it is 
impossible for all loops to have 4 A matrices running around them. This is simply the 
statement of the fact that it is impossible to put a completely flat lattice on the sphere. 
Positive curvature defects have to be introduced to close the surface. 

Since the order contribution is trivially zero, it is interesting to investigate the 
hrst non-trivial order: the order contribution. We will give some qualitative arguments 
that the flat case is described by N fermions with an equidistant spectrum and, from this 
starting point, calculate the number of flat graphs in the hrst 2 orders in 1/N. 

One may notice that, in the large N limit, all the factorials and products of differences 
of various highest weights cancel (assuming all 4 groups of h’s to be distributed in the same 
way). One can hypothesize that this will be true even for the next order in 1/V^ which 
describes the graphs with the topology of a torus (we will check this assumption below). 
We are then left with the partition function of fermionic type: 


Z{X) = A^(^-b/2 (6.4) 

hiA>h2-..>hN 
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A standard calculation for the free energy /(A) = l/A^^lnZ(A) gives: 


N OO 

/(A) = -J2 ln(l - A’") = -iV -2 ln(l - A’") + 0{N-'^) = + C»(iV-4), 

n=l n=l 

(6.5) 

where 77 (A) is the Dedekind function. The order N'^ contribution, which counts the number 
of flat graphs with spherical topology, is zero since, as already stated above, no such graphs 
exist, at least not without defects. The next order counts the number of flat graphs that 
can be htted onto a torus. It is easy to verify this calculation by directly counting the 
number of possible graphs. A general graph on the torus consists of ??r x n squares glued 
into a rectangle {m columns and n rows). Opposite sides are then glued together: hrst we 
glue together the two sides of the length n and then, with m possible twists, the two sides 
of length m. The symmetry factor is l/{mn). We obtain: 

00 00 

/(A)= y; —A’"" = - V ln(l - A"), (6.6) 

n,m=0 n=0 


which coincides exactly with our result obtained from the sum over highest weights treated 
as free fermions. Since completely flat graphs exist only for the torus topology we can be 
sure that there exist no higher order contributions. Let us note that this calculation on 


the torus is similar to what we would have for QCD2 on a toroidal target space [|13 


We conclude that, in the limit of flat graphs, the highest weights play the role of the 
energy levels of N free fermions. This fact may be useful for the understanding of the 
mechanism of the flattening phase transition. 

We are now in a position to discuss the real problem of interest. We have observed 
that, for models 


and (| 6 . 1 |) , where both the vertex and face coordination numbers 
are even, the simple saddle point equation is valid. In fact, all the simple models with 
even face and vertex coordination numbers which we have investigated (for brevity we do 
not discuss them here) can be solved correctly with the simple saddle point formulation. 
We thus hypothesize that if we restrict our attention to the “even-even” models (models 
where the coordination numbers for both faces and vertices are even), the destabilization 
of the saddle point discussed at the end of section 5 will be avoided and we can correctly 
interpolate all the way from the A = B = J case ( b. 8 |) to the flat lattice (|6.1|) using the 
saddle point approximation. Any results we obtain can be checked against perturbation 
theory to verify this assumption. For the interpolating model we thus choose 


A = B = 


Vb 

0 


0 

—a/6 


(6.7) 
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so that all odd traces equal zero, calculate the character from (|2.4|) (see Appendix) 


X{h}{A) = X{h}{B) = X{t^}iP)X{h!i^^{h) sgn[JJ(h^ - h°)], 

i,j 


(6.8) 


and arrive at the partition function 


{h‘^,h°} i 


Ajh’-f A{h‘f . A , |w(w+i)+^y^ ().;+>n 

HijiK-h’) ^N> 


X 



h° — 1 


(6.9) 


where — X{i^}(^)^(^)/(^(/5)^('^)) is an Itzykson Zuber integral, with b = 

and I{B^l^b) is dehned similarly. The partition function is a generalisation of (|5.8|) the 
difference being the two Itzykson Zuber integrals. Indeed, setting 6 = 1 we recover ( |5.8|) . 
For convenience we introduce the notation ig = -^Tr b^ so that ig = t 2 q = t^g- In this 
notation, flattening of the lattice corresponds to setting t 2 ^ 0 and tg = 0 for g 7 ^ 2. In 
fact, in complete analogy with the “rainbow” —> “cigar-like” transition in the gaussian 
model, it is only necessary to set tg = 0 for g > 2 to flatten the lattice. For convenience, 
we dehne F{h^) — 2-^I{^,b). It is then easy to derive the following pair of equations 
which in principle completely describe the model: 


qig = —Tr 6^ = - / 

N q J 2Tii 


i 


dh' AA = _ in(Aft) 


In 


h 


h-b 


2F{h). 


( 6 . 10 ) 


( 6 . 11 ) 


The hrst comes from a generalisation of (^) and the second is the saddle point equation. 
The whole complexity of the random to flat transition is succinctly encapsulated in these 
two equations. As in the toy, gaussian model, we can capture the transition by turning on 
the hrst three coupling constants tg. However, this system of equations, though simpler 
than the original statistical mechanics problem (0), is still highly non-trivial. As it is 
a kind of Riemann-Cauchy problem, the difficulty is as always to guess the analytical 
structure of the solution. 


7. Conclusions and discussion 

In this paper, we developed a character expansion technique for a new kind of ma¬ 
trix model describing dually weighted planar graphs. We hrst reduced the N'^ degrees of 
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freedom of the original matrix to the N highest weights specifying the irredncible repre¬ 
sentations of U{N). This allowed ns, with some precantions, to apply the saddle point 
method in the large N limit, redncing the problem to a set of integral eqnations. 

We showed how to solve these integral eqnations in a nnmber of known 1-matrix 
models and also fonnd a new resnlt: we are able to calcnlate the nnmber of graphs having 
only even nnmbers of neighbonrs for both original and dnal vertices. 

We also nnderstood the limiting case of the completely flat lattice. It is described by 
a system of free fermions (the highest weights) with an eqnidistant spectrnm. This allowed 
ns to reprodnce correctly the partition fnnction of regnlar graphs with toric topology (the 
only genns that can be realized from a completely flat graph). 

The most important physical qnestion still to be addressed is the description of the 
transition from completely random planar graphs (describing the 2D gravity) to the regnlar 
(flat, in onr terminology) lattice. 

We have not yet solved the corresponding integral eqnations. It is not an easy task as 
they are eqnivalent to a complicated Canchy-Riemann problem, and the solntion involves 
some non trivial gnesses abont the analytical strnctnre of the nnderlying fnnctions. 

Let ns specnlate possible physical pictnres for the flattening phase transition. There 
are three different scenarios to consider: 

I. The flattening conld take place for a hnite effective conpling constant in front of the 
(cnrvatnre sqnared) term in the 2D gravity action. This means that the characteristic 
flat size of an almost flat piece of a graph diverges at some hnite critical conpling. This 
wonld be the most interesting scenario as it wonld mean the discovery of a completely new 
nniversal critical phenomenon. 

An argnment against this pictnre is the fact that the conpling is dimensionfnll, 
thns containing inverse powers of the cnt-off. This is however a completely pertnrbative 
argnment since the same reasoning conld be applied to the 3D Ising model described as 
(f)'^ scalar held theory. Here, the interaction term is also dimensionfnll, bnt it nevertheless 
exhibits a non-pertnrbative phase transition. 

A more serions objection to a transition at hnite ehective conpling is the presence of 
macroscopic excitations on the backgronnd of a regnlar lattice generated by only a very 
small nnmber of cnrvatnre defects. For example, on the sqnare lattice, the introdnction 
of fonr vertices of coordination nnmber 6 and fonr of coordination nnmber 2 allows a 
baby nniverse of arbitrary size to grow ont from the hat lattice. This is in complete 
analogy with the tree-like strnctnres stndied in the toy ganssian model of section 4, where 
the introdnction of jnst one ^3 defect and one t\ defect is enongh to create a whole new 
branch. 
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2. The flattening could appear only in the limit of the infinite B? coupling. In the light 
of the previous arguments this is the most likely scenario. Another argument was given 
in using the methods of Liouville theory, though again the argument was completely 
perturbative and thus unsatisfactory. This hypothesis also seems to be in agreement with 
numerical simulations . Even if this scenario turns out to be correct, the asymptotic 
approach to the flat lattice could contain some interesting scaling behaviour and is worthy 
of study. 

3. The third scenario is an extended version of the second one: we could have a 

flattening phenomenon for the coupling of the order of course depend 

upon the type of regular lattice, and the critical exponents, if any, will be dependent on 
the symmetry of the lattice; triangular or square. In this case the phase transition would 
be better identified as a crystalization transition. The universal critical properties would 
then depend on the symmetry of the ” crystal”. 

One could think of an analogy with the Berezinski-Kosterlitz-Thouless phase tran¬ 
sition: the discrete curvatures on the lattice could be identified with quantized coulomb 
charges in two dimensions. This is indeed the picture in the conformal gauge of the con¬ 
tinuous 2D gravity action. 

Let us also make some technical comments. The character expansion method can be 
successfully applied to a more general type of matrix model, having the following action: 

S{M) = tr{W{M)+ V{AM)) (7.1) 

instead of just W{M) = as in the present paper. The character expansion coefficients 
are more tricky, as in |^ , but there are no real obstacles: we still deal with only N highest 
weights as the principal degrees of freedom. 

One can imagine the solution of a generalized 2-matrix model, such as 

Z = j exp (Tr(l/(L) + V{M) + W{LM ))), (7.2) 

with an arbitrary function W{x), not just a merely linear one, or even a similar multi¬ 
matrix chain. 

The most ambitious step would be to put ’’matter” on a random, but gradually flat¬ 
tening, lattice. For example, one can consider the matrix action: 

S{M) = tr{W{AL) + V{AM) + L'^ + + cLM). (7.3) 

It describes the Ising model on DWG. This model provides the interpolation between two 
solvable cases: the Onsager solution for a regular lattice and the Ising model on dynamical 
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random graphs |]K| . li W V it includes the non-zero magnetic held, which is still an 
unsolved case for the regular lattice. Unfortunately, the character expansion would be 
much more complicated object in this case as it would contain some non-trivial Klebsch- 
Gordan coefficients. 

In any case, we believe that the proposed approach could be fruitful for attacking 
many new combinatorial problems in 2D statistical mechanics and held theory. 


8. Appendix 

Below we present explicit formulea for some characters derived from the dehnitions in 
section 2. 


8.1. A = 1, the unit matrix 

In this case the character is just the dimension of the representation. The easiest way 
to derive this is to take the limit as e ^ 0 of the character formula, equation (|2.4| ) with 
cik — In this case 

where c is the numerical constant c = ■ 


8.2. Am, dehned by ^Tr(A^) = 4,^ 

The traces of all positive powers of Am are zero except We will sketch the 

derivation for the case m = 2. It is easy to generalize the derivation for arbitrary m. Using 
the second dehnition for the character, equations (|2.(j|) (O), we have 


Pk 


(fc/2)! 

0 




k even and non-negative; 
otherwise 


( 8 . 2 ) 


If we substitute this into the determinant, we obtain a matrix structure in which every 
other entry in a row is zero. By interchanging rows and columns, the determinant can be 
put into block diagonal form, with one block for the h^’s and the other for the h°’s. The 
powers of N/2 factor out and if we then factor out the product (ni(^)K“V~)0 
entries in the diagonal blocks become polynomials of ascending order in or h°, and the 
block determinants reduce to Vandermonde determinants. Taking into account all of the 
sign changes from reordering the rows and columns we obtain: 


X{h}{M) = c 



A(h^)A(h°) 


sgn[]q(/i' 



(8.3) 


20 










where c is a numerical constant. For general m we have: 


f N\ 

X(m(-4„)=c (-j 


m— 1 


n 





)! 


sgn[ Yl 

0<ei<e2<{m—l) i,j 


(8.4) 


In this case, the integers h factor into m groups of ^ integers with e = 0,1,..., (m — 1) 
denoting their congruence modulo m. 


8.3. Af,h defined by Tr = 0 for odd q 

Only the even powers of A^t are non-zero. The matrix A^b can be dehned by an ^ 
by ^ matrix, 6, (eigenvalues bk) as follows: 


Abb — 


b 

0 



(8.5) 


Again we just sketch the derivation. Calculating the determinant in equation (^]^) and 
rearranging the columns we notice that we can write it (up to some sign factors from the 
interchanges) as 


*<(»,■) (of) 


bh^ 


b^^ 0 

6 ^^ -b^” 

= (- 2 ) 2 

0 b^° 


(8.6) 


where, for notational convenience, we denote by b^ the ^ T whose elements 

are bfi. The Vandermonde in equation (|2.4| ) is just the special case of the above result, 
( F^) , with = 2j — 2 and h° = 21 — 1. Including the sign factors neglected earlier we 
arrive at the formula 


X{h}{Abb) = sgn[JJ(h^ - /i°)]. 


(8.7) 


To obtain the character for the ff matrix introduced in section 5, we set 6 = 1 and obtain 
the expression given in (|5.1| ). 

It is easy to generalize to higher order cases. For example, to study the case where 
only every third power of the matrix Abbb has a non-zero trace, we start with an ^ by ^ 
matrix, 6, and dehne 


Abbb 


boo 

0 ub 0 
0 0 Lu'^b 


( 8 . 8 ) 


where uo is the third root of unity. This time the character factors into three characters, 
one for each of the congruence classes, modulo three, of h. 
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